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Abstract 

We introduce a novel class of spatial complex variable function Gaussian (CVF-Gaussian) beam, 
which is the product of an arbitrary analytic complex variable function and a Gaussian function, 
in strongly nonlocal nonlinear media (SNNM). The CVF-Gaussian beam rotates generally during 
propagation. By choosing the input power of the beam, we can obtain the CVF-Gaussian breather 
or the CVF-Gaussian soliton. We reveal that stable CVF-Gaussian beam can exist in SNNM with 
different forms, including rotating dipole, and rotating elliptic donut. A distribution factor which 
is the parameter for the description of the transverse distribution of the CVF-Gaussian beam is 
discussed. 
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Recently increased attention in the study of self-trapped optical beams in nonlocal non- 
linear media, specially in strongly nonlocal nonlinear media (SNNM), is due to the ex- 
perimental observations of nonlocal spatial optical solitons (NSOSs) in liquid crystals [l| 

A 



n 



and lead glasses |2j, as well as many interesting theoretical predictions ^, 4, 5, (J 
variety of NSOSs, including dipole solitons [7], vortex solitons [8j, Hermite-Gaussian soli- 
tons (HGSs) and Laguerre-Gaussian solitons (LGSs) [o], Ince-Gaussian solitons (IGSs) 10], 



azimuthons [111 ] . ellipticons 



12| and rotating solitons Q have been presented. Ex- 



ed "walk-off" and steering of the NSOS [15] , 



perimenters have observed the voltage-control 

both the long-range [l(| and the short-range |l7l. Il8l] interactions between the NSOSs, the 



nonlocality-controlled interaction of the NSOSs 19j, as well as boundary forces of the media 
on the NSOS {20]]. Multipole solitons in SNNM, including dipole, tripole, quadrupole, and 
necklace- type solitons were also experimentally observed |21j ]. 

In this letter, we introduce a novel class of spatial self-trapped optical beam in the SNNM, 
the so-called spatial complex variable function Gaussian (CVF- Gaussian) beam which is 
the product of an arbitrary analytic complex variable function and a Gaussian function. 
Depending on the input power of the beam, the CVF-Gaussian beam can be either a soliton 
state or a breather state. 

Our analysis is based on the nonlocal nonlinear Schrodinger equation (NNLSE) for the 
amplitude $(r, z) of a paraxial laser beam propagation in a focusing cubic media with a 
transverse nonlocality [o, Q : 



id^/dz + ^A_l$ + kAn/n § = 0, 



(1) 



where z is the longitudinal coordinate, fi = I /(2k), k = ujtiq/c is the wave number in the 
media without nonlinearity, n is the linear refractive index of the media, Aj_ is the two 
dimensional transverse Laplacian operator, An = n 2 ff^° R(r — ri)|$(ri, z)\ 2 d 2 vi is the 
nonlinear perturbation of refraction index, n 2 is the nonlinear index coefficient, r and r x 
are the two dimensional transverse coordinate vectors, and R is the normalized symmetrical 
real nonlinear response function of the media. 



In the case of the strong nonlocality 



ity 



121 ] . if the response function R(r) is twice 



differentiate at r = 0, Eq. (TjQ) can be reduced into [J, [17|, [22| the Snyder-Mitchell model |6|] 



id^/dz + /iA ± ^ - n 2 /(2n )kjP r 2 ^ = 0, 



(2) 
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and the functions $ in Eq. ([I]) and \I/ in Eq. (|2j) are related by [17] $ = 
\I/exp [ifcr^-Ro-Fb-z/wo — in2/(2no)kjPo Jq Z / |^( r ij | 2 rf cZ 2 ri] , where r = |r| is the ra- 
dial coordinate, 7 is the material parameter associated with the response function R(r) {4], 
Pq is the input power at z = 0. The Gaussian beam solution of Eq. is ^/q — 
\/P^ey^)[i6(z)\/[^/Kw(z)]ex\){ — r 2 /[2w(z) 2 ] + zc(z)r 2 }, where w(z) is the beam width of 
the Gaussian beam, c(z) represents the phase-front curvature of the beam, 8(z) is the phase 
of the complex amplitude. They are given by respectively, w(z) = w (P c / P sin 2 f3 z + 
cos 2 (3 z) 1/2 } c(z) = k[3 w 2 (P c /P - 1) sin2/3 z/[4w 2 (z)], and 



9(z) = — arctan ( a/ P c / Pq tan (3 Q z) , (3) 

where wq is the initial beam width at z = 0, P c = uq/ ' (k 2 ^n2W ) is the critical power for the 
soliton propagation, (3 = ^/n 2 / yPo/n . To obtain CVF-Gaussian beam solution of Eq. (j2J), 
we suppose its solution has a form, 

* = # F (r,z)# G (r,z). (4) 

Substituting Eq. (J4j) into Eq. (j2J), one can obtain, 

.<9# F 



<9z 



2// 
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2ic(z) l x —£- + y— £. =0. (5) 



-w 2 (,z) J\ & cty 

In the transverse plane perpendicular to z, we introduce a rotating coordinate system 
(x',y',z), x' = wo[x cos(i?(z)) + y sin( , #(,?:))]/w(,2!), y' = Wo[y cos(i?(z)) — x sin( , i9(2;))]/w(z), 
where describes the rotation angle of the rotating coordinate, dd(z)/dz denotes the 
angular velocity (AV, the angle rotated per unit propagation distance). In the rotating 
coordinate system, Eq. (jSJ) can be written as, 



w 2 (z) dz 
2jj_( ,d* F ,dV F \ 

~ l^z)V^ + y W) = ^ (6) 

where A' ± = d 2 /dx' 2 + d 2 /dy' 2 , L' z = —i(x'd/dy' — y'd/dx') is the z component of the 



angular momentum operator, equivalent to one in quantum mechanics 23J. The solution of 
Eq. ([6j) is found to be 

*f(v',z) = f(^f), *(z) = T0{z), (7) 
3 



where /(•) denotes an arbitrary analytical function, b 7^ is an arbitrary real parameter 
whose physical meaning will be discussed later, 9(z) is given by Eq. (j3J), and " + (— )" 
presents the rotation direction which accords to the right-hand (left-hand) rule relative to 
the propagation direction. From Eqs. (j3j) and (j7|), one can obtain the solution of Eq. (j2j) in 
the laboratory frame: 



ty(r, z) = — — exp 



+ ic{z)r 2 + iO(z) , (8) 



w{z) \-2w(z) 2 

where Co is a constant which can be determined by ff™ |\&(r, 0)\ 2 dxdy = Po, (± = (x ± 
iy)/[bw(z)]exp[^pi$(z)]. The solution is called CVF-Gaussian beam due to the fact that 
/(C) is an arbitrary analytical complex variable function. The structure of the beam is 
determined by the product of the function /(C) and a Gaussian function. Solution (jSJ) 
represents an optical beam that rotates with an AV, d$(z)/dz, in the laboratory frame 
(x, y, z) and is stationary in the rotation frame (x', y', z). Observed in the laboratory frame, 
the trajectories of the corresponding points of the beam during propagation rotate generally. 

To discuss the physical meaning of the parameter b, we define the mean radius of an optical 
beam as [24] W(z) = 2 ((r 2 ) - (r) 2 ) 1/2 , where (X) = Jj^ X\^\ 2 dxdy/ Jj^ \^\ 2 dxdy. From 
the axis-symmetrical solution (jSj), it can be obtained (r) = 0, which means that the position 
of the CVF-Gaussian beam center is always located on the propagation axis (z-axis). By 
use of the Taylor expansion of /(C), then, the mean radius of the CVF-Gaussian beam is 
found to be 

2 _ Mzf ~ \A m \\m+l)\ 

p 2^ b 2m > ^ 

u m=0 

where A m = d m f(Q/d( m \^ = o/m\ (m = 0, 1,2, • • •) are the coefficients of the Taylor se- 
ries of /(C) at C = 0, the input power P = J^ =0 \ A m \ 2 m\/b 2m , and ml is the fac- 
torial of an integer m. W characterizes the distribution of the CVF-Gaussian beam 
from the beam center and provides a way to capture the beam's scale or degree of be- 
ing spread out. The derivative of W 2 with respect to b can be written as dW 2 /db = 
-Aw{z) 2 /{bP 2 ) [P J2™ =0 \A m \ 2 m 2 m\/b 2m - ( £~ =0 \A m \ 2 mm\/b 2m ) 2 } . By mathematical in- 
duction, we can prove that dW 2 /db = if /(C) = A or /(C) = Ai(, where A ^ and 
A\ 7^ (It means that the function /(C) is a constant or direct proportional to (.), otherwise 
dW 2 /db < 0. This means that the mean radius of the CVF-Gaussian beam is independent of 
b for the former case, and is a monotonically descent function of b for the latter case. There- 
fore, we can find that b, which is nominated a distribution factor, is the parameter for the 



description of the transverse distribution of the CVF-Gaussian beam, and conclude that the 
distribution is independent of b if the function /(£) is a constant or direct proportional to (, 
otherwise generally for another given /(C), the distribution is farther from the beam center 
for a smaller b and is clustered more closely around the beam center for a larger b. When 
b — > oo, we have \P — > f(0)^a [see, Eq. PJ]. \l/ will, then, equal the Gaussian function or the 
null function depending upon whether /(0) ^ or not. Thus we have the following limit: 
W b ^°°) w(z), if f(0) 7^ 0; or W fe ^°°> 0, if /(0) = 0. Figure [1] gives an example of the 
transverse normalized intensity distribution of the CVF-Gaussian beam when f(() = cos(£) 
or /(£) = Pu(C) for different b, where H n (Q is the nth-order Hermite polynomial. Such a 
cosine CVF-Gaussian beam or fourth-order Hermite CVF-Gaussian beam is clustered more 
closely around the beam center as b increases, and will approach to a Gaussian distribution 
when b — > oo since cos(£)|(=o = 1 and i?4(C) |c=o = 12. In fact, the cosine CVF-Gaussian 
beam becomes almost the Gaussian function when b = 10, as shown in the inset of Fig.[T](d). 

The CVF-Gaussian beam can exhibit the different properties for different input power 
Pq. Equation (J9j) shows that the mean radius of the CVF-Gaussian beam is proportional to 
the beam width of the linked Gaussian beam when the function f(() and the distribution 
factor b are given. The CVF-Gaussian beam, therefore, evolves in the similar way the 



Gaussian beam does [4]. When P < P c [the expression of P c was given under Eq. (j3J)], 
the pattern of the CVF-Gaussian beam initially broadens and then narrows, undergoing 
a periodic evolution; whereas the reverse happens and the pattern initially narrows when 
Po > Pc- These two cases are the CVF-Gaussian breather states. When P = P c , Eq. (jSJ) 
can be reduced as 



V a (r,z) = ^ ' exp ( -— + i(3 oZ ) , (10) 



Cof(e±>) ( r< 
ex P -TT-i 

where = (x ± iy)/(bw ) ex.p(±i(3 z). This is a CVF-Gaussian soliton, and diffraction 
is exactly balanced by nonlinearity in this case. Its transverse distribution preserves but 
rotates as the CVF-Gaussian soliton travels along the propagation axis. 

From Eqs. (E]) and ([3]), the AV of the CVF-Gaussian beam can be obtained 

_*&{?)_ , [P~ c w 2 /3o 

dz _± VPo^ 2 W 1 j 

Equation (TTT]) shows that Q is only dependent upon the material parameters and input 
power, and independent of the complex-variable function / and the distribution factor b, that 
is, independent of the concrete form of the CVF-Gaussian beam. Q is inversely proportional 



to w 2 , the squared beam width of the linked Gaussian beam, so fl behaves periodically as 
w does 4j. If Pq < P c , Q decreases initially, vibrating between a minimum Pq/(P c zr) and 
a maximum 1/zr] whereas if Pq > P c , the reverse happens, and Q increases first, vibrating 
between a maximum P /(P c zr) and a minimum 1/zr, where zr = kw 2 , is the Rayleigh range. 
When Pq = P c , Q of the CVF-Gaussian soliton is a constant 1/zr. 

The analytical prediction above is tested by comparison with the numerical solution 
of Eq. ([I]). The material response is supposed to be the Gaussian function 
R{r) = l/(27ru£)exp [-r 2 /(2w 2 J] , where w m is the characteristic length of the material 
response function, a = wo/w m denotes the degree of nonlocality for the beam in nonlocal 
nonlinear media |4J. The less a is, the stronger the nonlocality is. The propagation dy- 
namics of the cosine CVF-Gaussian soliton, and the fourth-order Hermite CVF-Gaussian 
soliton are presented in Fig. [2], respectively. Figure [2] demonstrates excellent agreement 
between the exact analytical CVF-Gaussian solution of the Snyder-Mitchell model and the 
numerical simulation of Eq. (CQ) with the corresponding inputs given in the caption of Fig. [H 
Figure [2(a)- (h) shows that the beam pattern of the cosine CVF-Gaussian soliton is com- 
posed of rotating dipole-mode soliton. For b = 0.5, the pattern of the fourth-order Hermite 
CVF-Gaussian soliton consists of a bright elliptic ring with a dark elliptic spot, thus it seems 
appropriate to refer to this kind of ring soliton as "elliptic donut mode" [251 ] . 

A better understanding of this phenomenon of the visible rotation for the optical beam 
would be provided if an analogy between the rotating beam and the mechanically rotating 
object |26[] was drawn. Being similar to mechanical attributes, the optical beam possesses 
the AV, the moment of inertia that is proportional to the total beam equivalent mass, 
which equals the beam power, and to the squared mean radius, and the angular momentum 
(AM) that is the product of the AV and the moment of inertia. Just like the movement of 
the rigid body in mechanics, this shows that if the beam is narrowed or broadened while 
keeping the total "mass" , its AV increases or decreases. In this way, the result of Eq. (fTTjl 
could be clarified physically. From the definition of AM [141 ] . A = Im Jj^ $*(r x V$)<i 2 r, 
one can obtain the AM of the CVF-Gaussian beam is ^^ = s{\A m \ 2 mm\ /b 2m ) . The inverse 
proportion between the AM and distribution factor b could also be explained qualitatively 
via this analogy. Since the mean radius of the CVF-Gaussian beam is inverse proportional 
to b and the AV is independent of b, as discussed above, it is easily understood why the AM 
decreases while b increases with the total "mass" preserving. 
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Although the visible rotation comes from the non-zero AM, not all the beams with non- 
zero AM will rotate. It can be observed that the vortex LGS and the vortex IGS 10], 



which are also the solutions of the SNNM d2J), do not rotate during propagation, but both of 
their AMs are non-zero. The physics behind the rotation phenomenon is the inhomogeneity 
of the energy flow within the beam cross section. Extending the linear analysis 27j to the 
nonlinear case, we can have d\§\ 2 /dz = — Vj_ ■ Sj_, which follows directly from the Eq. (pQ), 
where Sj_ = i[$*V_i_$ — $Vj_$*]/(2fc) is the transverse Poynting vector component that 
represents the energy flow within the beam cross section, Vj_ = (d/dx, d/dy), and * denotes 
the complex conjugate. According to the definition of the divergence of a vector field, we have 
d\$>\ 2 /dz = — lim §„ Sj_ • dO/V, where the quantity Sj_ • dO is the scalar flux of the field 
through the closed surface 5 that surrounds a point (x, y, z) and encloses a domain D with 
volume V . When the energy flow is homogeneous, the energy outflow will equal to energy 
inflow through the closed surface, then <f„ S± ■ dO = and the intensity pattern is constant 
during the propagation. If the energy outflow is less than the energy inflow, §„ Sj_ • dO < 0, 
and d\Q\ 2 /dz > 0, hence the intensity grows in a point, otherwise the intensity falls down 
if the energy outflow is more than the energy inflow. It can be found that Sj_ 7^ but 
Vj_ • = for the vortex LGS and the vortex IGS, representing their energy flow within the 
beam cross section is homogeneous, therefore they do not rotate during propagation. For the 
CVF-Gaussian soliton, however, we have V ± ■ S± = 4ilm[f((¥)(¥*df(c¥*)/d(¥*]/(bw 2 ). 
Such inhomogeneity of the energy flow within the beam cross section results in a transverse 
energy redistribution over the beam cross section, and as a result the visible rotation of the 
CVF-Gaussian soliton can be observable. 

One more important issue should be talked about in the end. Numerical simulations 
of Eq. ([1]) when the initial beam is perturbed by noises 28] have proved the stability of 
the solution (JSJ). The detail of the simulation will not discussed here because of the length 
limitation. 

In conclusion, we have provided analytical and numerical evidence for the existence of a 
novel class of the CVF-Gaussian beams, the CVF-Gaussian breather and the CVF-Gaussian 
soliton, stabilized by strongly nonlocal nonlinearity. The CVF-Gaussian beams will rotate 
during propagation, and the angular velocity is generally a function of the material param- 
eters and the input power, independent of the concrete form of the CVF-Gaussian beams. 
The distributions of the CVF-Gaussian beams are generally dependent upon the distribution 
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factor 6, and the distribution is farther from the beam center for a smaller b and is clustered 
more closely around the beam center for a larger b. The stable CVF-Gaussian beam can ex- 
ist in SNNM with different forms, including rotating dipole, and rotating elliptic donut. The 
physics behind the rotation phenomenon is the inhomogeneity of the energy flow within the 
beam cross section, which results in the energy redistribution over the beam cross section. 
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FIG. 1: (Color online)The transverse normalized intensity distributions at z = plane 
for the cosine CVF-Gaussian beam ^(r, 0) = \J Pq sec h(l /b 2 ) / (y/irwo) cos(£)exp [ — 
r 2 /(2wQ)] [(a)-(d)], and the fourth-order Hermite CVF-Gaussian beam ^(r,0) = 



y/Pob A /[4:Wo ^3^(128 + 966 4 + 3b 8 )]H A (() exp [ - r 2 /{2wl)] [(e)-(h)]. The different columns 
represent the different b, as given in the top. The inset in (d) is the transverse normalized intensity 
distribution of cosine CVF-Gaussian beam (solid line), and that of the Gaussian beam (dotted 
line) at y = 0. 
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FIG. 2: (Color online) Propagation dynamics of the cosine CVF-Gaussian soliton [(a)-(h)] and the 
fourth-order Hermite CVF-Gaussian soliton [(i)-(p)] in the Gaussian-shaped response material. 
(a)-(d) and (i)-(l): transverse normalized intensity distributions from the analytical solution, (e)- 
(h) and (m)-(p): distributions from the numerical simulation (a = 0.02). The different columns 
represent the different propagation distances, as given in the top. The parameters are chosen as 
P /P c = l,b = 0.2 [(a)-(h)] and b = 0.5 [(i)-(p)]. 
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